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We theoretically investigate interactions between non-degenerate mechanical oscillators mediated 
by a time-dependent cavity field. We obtain a reduced master equation valid for all optomechanical 
systems operating in the weak coupling regime. This master equation includes all forms of deco¬ 
herence and back-action due to the dissipation of the field mediating the interaction. We apply 
the master equation to study two resonant coupling schemes within a rotating-wave approximation: 
the beam splitter Hamiltonian and the two-mode parametric amplifier. In both cases, the effective 
unitary interaction can be made arbitrarily strong compared to the decoherence due to dissipation 
of the mediating field by choosing appropriate detunings. 
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I. INTRODUCTION 

The ability to control a physical system at the quan¬ 
tum level has, whenever it was achieved, profoundly influ¬ 
enced physics and many of its applications. Two promi¬ 
nent examples, the development of lasers and their appli¬ 
cation to cool atomic gases, have allowed to use light as a 
tool to perform a wide range of tasks and lead to improve¬ 
ments of precision measurements by orders of magnitude. 
In the last few years, successes in cavity optomechanics 
have lead to similar control over motional degrees of free¬ 
dom of mesoscopic and microscopic oscillators [1]. 

Experiments have successfully demonstrated quan¬ 
tum effects involving a single mechanical oscillator EHz] 
and theoretical studies point out the suitability of op¬ 
tomechanical systems to observe truly macroscopic non- 
classical behavior [HHia- While the vast majority of the 
experimental results involve a single mechanical mode, 
typical solid-state resonators exhibit a multitude of res¬ 
onances over a wide range of frequencies [13 El]. In sys¬ 
tems where the mechanical modes consist of the collective 
motion of a trapped atomic cloud giisi one can engineer 
non-degenerate optomechanical modes by trapping sev¬ 
eral clouds in an optical superlattice m- Numerous non¬ 
degenerate mechanical modes are thus a common feature 
of optomechanical systems, which feature, however, re¬ 
mains largely unexploited to date as most experimental 
results [HEIlEaEZ] involve nearly degenerate mechan¬ 
ical modes. 

While multimode optomechanics has been studied the¬ 
oretically for specific applications [IZIIMES], there has 
been no general formalism to analyze the possibilities of 
tuning non-degenerate modes into resonance. The aim 
of the present work is to provide that formalism, includ¬ 
ing the effects of quantum noise and measurement back- 
action of the mediating field. The resulting equations are 
applicable over the whole range of optomechanical sys¬ 
tems, from large scale interferometers for gravitational 


wave detection [28] to oscillations of ultracold atomic 
clouds |4] or optomechanical microcavities [U |29] . 

The effective interaction between non-degenerate me¬ 
chanical modes is suppressed by their frequency differ¬ 
ence. Together with the typically small optomechanical 
couplings this prohibits effective multimode physics in¬ 
volving modes with vastly different frequencies. However, 
other fields of physics successfully use auxiliary modes 
to resonantly couple non-degenerate states. In atomic 
physics, for instance, two laser beams appropriately de¬ 
tuned from a third state can drive a two-photon reso¬ 
nance between two non-degenerate atomic states [30] . In 
non-linear optics, two optical fields can be tuned into 
resonance by an appropriate pump creating a standing 
wave in a non-linear crystal m- The optomechanical 
analogue of these techniques is an oscillating intensity in 
a cavity field coupled to two non-degenerate mechanical 
modes and is the subject of this Article. 

A time dependent intracavity intensity is created by 
any anharmonic pump, the simplest example being the 
field created by two harmonic drives with distinct fre¬ 
quencies. The resulting cavity field can mediate a reso¬ 
nant coupling irrespective of the frequency difference of 
the mechanical modes. Recent theoretical studies have 
applied specific examples of this very effect for special¬ 
ized applications, such as the cooling of a low-lying me¬ 
chanical mode [25], mechanical phase-conjugation [24] or 
preparation of non-classical states m- The price for 
resonant mediated interaction between non-degenerate 
modes is additional decoherence due to the dissipative 
nature of the mediating field. The cavity field that pro¬ 
vides mechanical coupling also carries information about 
the mechanical state out of the cavity and this loss of 
unitarity is necessarily accompanied by noise entering the 
mechanical system [23] , 

In this Article, we derive a general and complete mas¬ 
ter equation for the behavior of two non-degenerate me¬ 
chanical modes interacting via weak coupling to a shared 
cavity field. Our treatment is complete in the sense that 
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all sources of decoherence, back-action and noise that are 
a consequence of the mediating field’s dissipation are in¬ 
cluded. Taking these effects into account is crucial for 
the operation of any optomechanical device at the quan¬ 
tum level and a prerequisite to the observation of non- 
classical effects. The resulting equation is general in the 
sense that we do not assume any specific detuning or 
rely on the system’s being in any particular parameter 
regime, apart from the weak coupling assumption which 
is accessible to virtually all optomechanical systems to 
date [1]. The derivation of our master equation does not 
rely on the mediating field’s being the mode of an electro¬ 
magnetic resonator. With a few simple modifications it 
will also describe interactions mediated by bulk-phononic 
modes [18] . Similarly it is straightforward to add individ¬ 
ual, Markovian heat baths for the mechanical oscillators 
and while we have neglected them for the discussion of 
the individual couplings, their effect is discussed in the 
final section. 

After deriving the general master equation, we will 
study the two resonant effective interactions described 
by the beam splitter Hamiltonian 

Ul = JBs(b\b 2 + hbl) (1) 

and the two-mode parametric amplifier 

= JpA (bih + , ( 2 ) 

where bj denotes the annihilation operator for phonons 
in oscillator j. The former allows the mapping of a quan¬ 
tum state from one oscillator to the other while the latter 
leads to non-classical correlations between the two oscil¬ 
lators, akin to parametric down-conversion in quantum 
optics. Most results obtained in single-mode optome- 
chanics consist of these two Hamiltonians resonantly cou¬ 
pling the cavity-mode with the mechanical mode via the 
matching of mechanical side bands with the cavity res¬ 
onance. Optomechanical cooling, for instance, is a real¬ 
ization of the beam splitter Hamiltonian together with 
the large dissipation of the electromagnetic field modes, 
while parametric driving is a realization of the para¬ 
metric amplifier involving mechanical and electromag¬ 
netic modes [7]- Realizing these couplings between two 
mechanical modes opens the door to more customized 
and precise sensing of forces and fields, the application 
of quantum computation algorithms to mechanical sys¬ 
tems m and the study of macroscopic decoherence [32] . 

Of crucial importance is the ratio between the strength 
of the mediated unitary interactions and the decoherence 
rates due to the mediating field being emitted from the 
cavity. This quantity is a measure for the ability of the 
effective interaction to create or maintain non-classical 
correlations between the two modes and thus a gauge for 
the classicality of the interactions [2311331 [34]. We find 
that the relationship between unitary coupling strengths 
and total decoherence rates depends on the system pa¬ 
rameters and the detunings of the pumps with respect 


to the cavity resonance, such that different choices of de¬ 
tunings turn out to be ideal for a given application. Most 
importantly, we find that even for systems where the cav¬ 
ity emission rate of the mediating field exceeds all of the 
involved mechanical frequencies there are regimes where 
the engineered interaction allows non-classical motional 
states or the transfer of nonclassical states between os¬ 
cillators. 


II. MASTER EQUATION 
A. System Hamiltonian 

Consider an optomechanical system consisting of two 
mechanical modes with natural frequencies uji and UJ 2 
coupled to the field of an electromagnetic resonator with 
single photon/single phonon coupling strengths gi and 
g 2 respectively. A sketch of a possible setup involving spa¬ 
tially separate membranes coupled to the field of an opti¬ 
cal Fabry-Perot resonator is given in Fig.[2 The mechan¬ 
ical modes need not be located on distinct oscillators and 
the mediating field does not have to be optical. The sys¬ 
tem could also consist of two non-degenerate drumhead 
modes of a capacitor in a superconducting LC circuit [7] 
with photons in the microwave regime or two samples of 
ultracold atoms trapped in an optical cavity [T6] . 

Neglecting all other modes of the mechanical system 
and the cavity, the Hamiltonian in the lab frame is given 


by 


= l-Lopt + 'Hm + Hi + 

(3) 

with 


Hopt = uJc^^a + r]{t)a^ + 

(4) 
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Hm = idj b\bj , 

(5) 


<^.c^ 

+ 

II 

(6) 


i=i 


and accounts for the intensity of the intracavity field 
being emitted to a reservoir at zero temperature at a rate 
K m- We have set h = 1, dropped constant terms and 
defined the input coherent field driving the cavity with 
the rate rj. The coupling of each oscillator to its thermal 
reservoir is neglected in favor of a focus on decoherence 
due to the dissipative nature of the field that mediates 
the interaction. Thermal effects are explicitly discussed 
later. 

To tune a desired coupling into resonance, we need the 
cavity intensity to oscillate in time. This is achieved by 
driving the cavity with a bi-chromatic pump, 

( 7 ) 
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FIG. 1: (Color online) (a) Sketch of a possible experimen¬ 
tal realization involving separate optomechanical membranes. 
The mechanical modes could also be on the same resonator, 
(b) A scattering process contributing to the coherent exchange 
of phonons in the resonant beam splitter interaction. The 
three other contributing processes result from time reversal 
and pump photons with frequency ujli scattering with oscil¬ 
lator 2. (c) A scattering process contributing to the paramet¬ 
ric amplihcation of the mechanical modes. The three other 
processes are obtained in the same way as in panel (b). 


The beat-note between the two drives will provide the 
necessary energy to tune either Hamiltonian 0 or 0 
into resonance. The density operator of the system p 
satisfies the master equation 


= -iiy^p] + KC{a)p, 
with the Lindblad superoperator 


C{a)p = apd) — -a^ap - 


L /s I ^ 

-pa' a. 


(8) 

(9) 


Without loss of generality we will assume uji > uj 2 and 
ujl 2 > ^Li‘ The next step is to unitarily displace the cav¬ 
ity field operator by a time-dependent complex function: 
d ^ Qf(t) + d. If we choose 

a{t) = (10) 


with 


( 11 ) 


where Aj = all source terms for d cancel and 

it only carries quantum fluctuations of the cavity field. 
Without loss of generality, we may choose the phases of 
Pj such that both aj are real. The displacement oper¬ 
ation introduces the term gj(bj + bj) driving 

the mechanical degrees of freedom. The time indepen¬ 
dent part of this term is a constant force that will shift 
the equilibrium position of each oscillator. Additionally, 
the beat-note between the two pumps will drive each os¬ 
cillator off-resonantly. The constant force does not affect 
the dynamics of the system. For the small mechanical 
linewidths typical of optomechanics [T] the off-resonant 
terms drive oscillator j into a coherent state with phonon 
number of the order of i—^ ^ weak couplings 

\0JLl-0JL2-0Jj\ ^ ® 

Qj his small oscillation can be neglected for all effects that 
are considered in the remainder of this paper. 

Similarly, we can also neglect the interactions 
quadratic in d and thus find in a co-rotating frame for 
the optical operators the Hamiltonian 

2 

'H = 'Hm+ T] 9jOik (6^ +S]). (12) 

j,k=l 


B. Tracing out the cavity field 


Tracing out the electromagnetic system will yield an ef¬ 
fective master equation coupling the mechanical modes. 
To perform the partial trace over electromagnetic degrees 
of freedom we restrict the optical Hilbert space to the op¬ 
tical vacuum and first-order coherences. This is formally 
a projection of the optical Hilbert space unto a lower¬ 
dimensional subspace as described below. Introducing 
Fock states of the optical subsystem |n)o and operators 

Pmn ~ o('^|p|'^)o5 
the reduction is formalized as 


Pmn = 0 m + n > 1. (14) 

By virtue of the displacement operation this is a reason¬ 
able approximation. The reduced density operator of the 
mechanical system becomes 


Pm=o(0|/)|0)o. (15) 


Note, however, that this identification can only be made 
after taking the trace over optical degrees of freedom in 
Eq. 0, as there are non-vanishing contributions from 
the one-photon sector. Taking this into consideration, 
we find the equation of motion 


dt ~ 


2 

+ E + Pj,pio] + /i.c.) , 

j,k=l 

( 16 ) 
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with h.c. denoting the Hermitian conjugate. The evolu¬ 
tion of the operator pio is governed by 


— i[?^ 77 T,, Pio] 'i ^ ^ ^ cyPlO- 


dpiQ 

dt 

7 ,/C=i 

_(17) 

Assuming that pm changes slowly on a timescale k ^ we 
can solve this equation by elementary transformations. 
The result is 


Pio = -i T] 


j,k=l 


P- 




Hi/2 + i(A/c + LUj) Hi/2 + i{^k ~ ^j) 


Pm- 

(18) 


This expression can be plugged into Eq. (16) to give the 


Here are 4x4 matrices 


M, = 


Di O: 

Ot D 


+ ’ 


reduced master equation for the mechanical system. In¬ 
troducing the vectors b = ( 61 , 62 , 61 , 62 ), where ^ de¬ 
notes the transpose, we obtain 

= -i[77m,/3m]+b'^Mibpm + /Omb'^M2b + b'''M3/)mb. 

(19) 

( 20 ) 

( 21 a) 
( 21 b) 
( 21 c) 
( 21 d) 


with 


rn±l -- 

L 1 J „„ 9mgn 2^^ ^^2 + i{Ai± OJn) 


Pt]mn ~ 3m9n E 


-D? 


] mn ~ 9mgn E 


T ^n) 


/d /2 i{^i T ^m) 


] mn ~9m9n E CtkCtl X 
k,l=l 

^i{Ai—Ak)t 




^K/2 + i{Ai±UJm) k/2 - i{Ai^ UJn) , 

( 21 e) 

2 

] mn —9mgn E CtkC^l X 
k,l=l 

^i{Ai—Ak)t ^i{Ak—Ai)t 


k/2 + i(A/ ± Um) - i(A/ ± cj^), 

( 21 f) 


and m^n e { 1 , 2 }. 

This is the general form of the master equation in the 
lab frame, including all off-resonant terms. It is easily re¬ 
duced to the system of a single harmonic oscillator (with 
uji = UJ 2 and gi = ^2 or 92 = 0 ) or a single harmonic 
pump (with Ai = A 2 and o^i = 0^2 or 0^2 = 0). As 
sketched in Fig. (b) and (c), mediated interactions in¬ 
volve at least two photon-phonon scattering processes. 
Thus all effective amplitudes appearing in Eqs. ( 21 ) are 
of second order in the single phonon/single photon cou¬ 
pling strength. 

The unitary evolution described by Eq. ([T^ is con¬ 
ceptually simple. Each pump shifts the resonance of the 
oscillators by the frequency 


lin([0('‘]ji - [O 2 ]u) = 6Q.i+Ri{t), 


( 22 ) 


where we have separated the time-independent spring 
shifts 6fti from time-dependent contributions easily 
found from Eqs. ( 21 ). The explicit form of the spring 
shifts is 


j;0 = „2 - ^i) _ alAk+^i) \ 

* ^'^UV4+(Afe-02^)2 K2/4 + (Afe+a;0V' 

(23) 


The mechanical resonances oscillate with small am¬ 
plitudes \Ri{t)\ around their means uJi + 6fti. Due 
to the weakness of optomechanical couplings we have 
\Ri{t)\^SQi uJi and the oscillation of their reso¬ 
nance frequency will not affect the dynamics of the sys¬ 
tem for sufficiently non-degenerate modes, quantified by 
|Ai — A 2 I ^ Additionally, the light-held also es¬ 

tablishes an optical spring between the oscillators, with 
spring constant given by 

ki 2 = -^-Im ([0+]i2 - [O^hi) , (24) 


where 2 (HO,i is the zero-point displacement of oscillator 
i. The re-scaling by the two oscillator lengths serves to 
give ki 2 the units of force per length but does not appear 
in the equations of motion, where this length has been 
absorbed into the dehnition of the operators 6 ^. This 
spring too consists of constant and time-dependent con¬ 
tributions. In contrast to the single-mode optical springs 
given in Eq. (22), it is now the time-dependent contri¬ 
butions that can create resonant couplings between the 
non-degenerate oscillators. 


III. RESONANT COUPLINGS 

In the following section we will focus on the cases where 
the time-dependence of the optical spring between the 
oscillator resonates with eigenmodes of the two-oscillator 
system. We will focus on a regime where the frequency 
difference between the two oscillators is large compared 
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to the effective coupling strengths, 


91920 ^ 10^2 

K 


< \uJl - UJ2\ 


(25) 


For such systems, a monochromatic pump does not lead 
to multimode behavior and the system can be treated as 
two separate single-mode optomechanical oscillators. In 
our case of two harmonic pumps, the beat-note between 
the two drives can tune a desired interaction into reso¬ 


nance. For such cases condition (25) assures off-resonant 


contributions to be suppressed such that they can be ne¬ 
glected within a rotating-wave approximation. 

The beat-note between the two optical pumps has to 
hit the resonances of the optically shifted mechanical fre¬ 
quencies uji + 6fti. This turns the resonance conditions 
into non-linear equations with respect to the optical de¬ 
tunings. The correction of the modified resonance con¬ 
ditions compared to their bare values are small, but we 
have computed them for completeness and they are given 
in the appendix. From here on, we assume that these cor¬ 
rections are taken into account and take them into the 
definitions of the mechanical frequencies, i.e. 


iUi —y uji -j- 


(26) 


This substitution has to be carried through in Eqs. ( [^ 
to be self consistent and does not lead to qualitatively 
different dynamics in the weak-coupling regime. 

Resonant coupling arises when the modulation fre¬ 
quency of the spring connecting the two oscillators 
matches the difference of the two individual oscillator 
resonances. 


Ai - A2 = (27) 

where we have introduced Sou = ooi — uj 2 for the differ¬ 
ence of the mechanical frequencies. The resulting inter¬ 
action - governed by Hamiltonian ^ - can be used to ex¬ 
change the quantum states of the two mechanical modes 
irrespective of their frequency difference. A diagram for 
one particular set of scattering processes contributing to 
this interaction is given in Fig. (b). The two clas¬ 
sical pumps driving the cavity establish the conditions 
for resonant exchange of phonons between the mechan¬ 
ical modes via scattering of intracavity photons. Alter¬ 
natively, the modulation of the coupling spring ki 2 can 
match the sum of the oscillator frequencies. 


Ai - A 2 = 2n (28) 

with O = being the average frequency of the bare 

mechanical modes. In this case the resulting interaction 
is a parametric drive of a collective mode of the two- 
oscillator system described by Hamiltonian (©■ An ex¬ 
ample for a process contributing coherently to this effec¬ 
tive interaction is given in Fig. (c). 

In the following, we will study the two resonant cases 
in detail. To this end, it is useful to transform into a 


(a) 

noj) 

A 


JL 


i 




‘1-0)2 




JL 

0)i^l+0)2 


i 




‘^0)1 


0 ) 


(b) 

P(^) 

A 


-- Zi -- ^ ^ 

0)i^2~^\ i^L2~^2 ^L2 ^^L2+^l 


FIG. 2: (Color online) Sketch of the intracavity photon power 
spectral density neglecting the cavity profile for clarity, (a) 
and (b) show the contributions of the two pumps in red (solid) 
and blue (dashed) respectively. The difference between the 
pump frequencies is arra nged in the beam splitter configura¬ 
tion specified by Eq. (27) with the pumps indicated by arrows 
representing delta-peaks at and (jUl 2 - Each pump creates 
a pair of sidebands for each mechanical oscillator, as indicated 
in the sketch. The overlapping sidebands are emphasized by 
the dashed lines and formally captured by Eqs. (33). The ex¬ 
act behavior of this quantity including the cavity profile can 
be obtained from Eq. (18) 


co-rotating frame for the mechanical operators, achieved 
by applying the unitary 




U = e 


(29) 


to the master equation (19). The transformed Hamilto¬ 


nian does not contain the free mechanical Hamiltonians 
Hm, but the frequencies now modulate the off-diagonal 
elements of the matrices appearing in the master equa¬ 
tion. 


[DfUn ^ (30) 

[OfUn ^ (31) 


A. Beam Splitter 

In this subsection we will study the setting where the 
effective interaction is given by the beam splitter Hamil¬ 
tonian 0 . This interaction leads to a coherent exchange 
of quantum states between the two mechanical modes at 
a rate Jbs- Applications include the cooling of a low- 
frequency mode by means of an auxiliary high-frequency 
mode, the study of mediated interactions [23] or quan¬ 
tum state mapping. In the time-domain the resonant 
coupling corresponds to a modulation of the strength of 
a spring connecting the two oscillators. It is instructive 
to discuss this condition in the frequency domain. The 
optomechanical interaction leads to side bands on the 
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cavity field, identifiable in Eq. (18) and associated with 
photon-phonon scattering processes. Each pump creates 
two pairs of side bands, one for each oscillator. Their 
optical frequencies are 


Ak i iOn 


(32) 


The resonance condition (27) implies 


= A2 — UJ2 — ^ ^ 

Ai + W2 = A 2 + Wi = A + fi, 

where we have introduced the central detuning 

^ ^ Ai + A 2 


(33a) 

(33b) 


(34) 


These identities represent overlapping side bands of dif¬ 
ferent oscillators, as sketched in Eig. When a photon 
coming from one pump mode leaves the cavity in the 
other pump mode, it will have interacted with both os¬ 
cillators and mediated the coherent transfer of a phonon 
between oscillators. The distinct overlapping side bands 
in Eqs. (33) provide two distinguishable pathways for 


phonon exchange and the coherent effective coupling rate 
is the sum of the amplitudes for each pathway. Addition¬ 
ally, there are four unmatched side bands with frequen¬ 
cies 


Ai T , (35a) 

A 2 -CJ 1 , (35b) 

A 2 +CJ 2 , (35c) 

Ai — co’ 2 . (35d) 


These unmatched side bands contain photons which have 
scattered phonons off a single mechanical mode. They 
carry information about the oscillators out of the cavity 
and are sources of dissipation and back-action. 

The resonant coup ling terms are found in the matrices 
using identities (p^ we find 


[of]i2 = [od2i = [02d:. = [02d;i = 


-gig2C^ia2 
f +i(AT^)’ 
(36) 


where the star denotes complex conjugation. The co¬ 
herent processes are contained in the imaginary parts of 
these coupling terms. 


JBs=Im([0+],2+[Or]i2)’ 


(37) 


where the two terms correspond to amplitudes for the 
scattering being mediated by the matched side bands on 
the red and blue side of the pumps respectively. The 
real parts of Eq. (36) together with resonant terms in O 3 


describe the dissipative dynamics. We find 
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ro±l = 9i92aia2K. 
y 3 J 21 - K2/4 + (A±ii)2' 


(38) 


Explicitly separating unitary and dissipative terms yields 
after some algebra the master equation 

^ = - z [77eff, Pr^] + r ((n + l)C{B)p + nC{B^)p) 

2 

+ E ((^J + + njC{b])p^ (39) 

i=i 

with the effective Hamiltonian 

'Heff = >^BS {b\b2 + & 1 & 2 ) ) (40) 

and the collective mode 



where we have assumed ai = a 2 = a and introduced 
Gj = gja for notational brevity. The chosen normaliza¬ 
tion of the collective mode yields non-canonical commu¬ 
tation relations but simplifies the coefficients appearing 
in the master equation. They are 


Jbs =Im 

ri = 


f?if? 2 (/z -f 2 zA) 

KV4 + fi2-A2-hmA, 

4GjKA{n + 5aj) 


(42a) 


To = 


r = 


ni = 


712 = 


n =- 


(k2/4 + (A - - (5w)2)(k2/4 + (A -k -k 5u})^) ’ 

(42b) 

4GlKA{n - Slu) 

-h (A -f- H — 6lo)‘^){k^/ i -I- (A — f) -|- 6u)‘^) ’ 

(42c) 

4 (ir (ir 2 A fl 


(k2/4 + (A - H)2)(k 2/4 + (A + fi)2) ’ 
k^/ 4 -h (A — H — (5a;)2 
4A(fi -f 5lo) ’ 

K?/A -f (A — fi -f (5a;)2 
4A(H - 5w) ’ 

k'^/A + (A - H )2 
4Afi 


(42d) 

(42e) 

(42f) 

(42g) 


The master equation (39) describes two resonantly cou¬ 
pled bosonic modes connected to three independent heat 
baths. The two heat baths described by coupling rates 
Vj and effective thermal occupations fij affect each os¬ 
cillator separately, while the bath described by T with 
effective occupation n acts on a collective mode of the 
system. 


1. Unitary Evolution 

Eor given side band resolution k/CI and optomechan¬ 
ical couplings, the maximum exchange rate between 
the two states Jmax is reached for A = ±{Q^hi/2). 
Eor typical experimental parameters the coupling is 
around 0.01 Cl [T] and can be tuned via the pump 
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^max 



FIG. 3: (Color online) Normalized coupling strength as a 
function of central detuning and side band resolution. The 
line gives parameters resulting in zero coupling, Jbs = 0. 
Other parameters are Gi = G 2 = 0.01 Q. 


strengths. The coupling strength as a function of central 
detuning changes its behavior depending on the side band 
resolution as plotted in Fig. For n <2 Cl there are two 
non-zero values of A resulting in no effective coupling, 
indicated by a line in the figures. This cancellation is a 
temporal analogue of a double-slit interference. The pro¬ 
cesses of phonon exchange via either matched side band 
(see Fig. and Eqs. (33)) can coherently cancel. This 
is, however, only possible if the photon coherence time 
exceeds the inverse separation of the two processes 

(cJi + • 



FIG. 4: (Color online) Scaled total decoherence rate F/Fmax 
as a function of detuning and optical dissipation rate. Other 
parameters are the same as in Figs. 


cavity resonance (A = 0) there is no effective spring cou¬ 
pling the two oscillators, Jps = 0. The light field only 
measures the mechanical system and its dissipation is 
pure measurement back-action of the light field on the 
mechanical oscillators m- Cavities with a dissipation 
rate smaller than the frequency difference of the oscilla¬ 
tors can resolve each heat bath independently, as plotted 
in Fig. 1^ In such systems it is possible to cool or am¬ 
plify each mode or the collective mode selectively and it 
is sufficient to consider only the dominant heat bath. 


3. Limits to Coherence 


The detuning resulting in the strongest coupling 
strength Jps is not necessarily the ideal choice to cou¬ 
ple two mechanical modes. For coherent manipulations 
on a quantum level it is necessary to minimize the influ¬ 
ence of decoherence. The total decoherence rate of the 
system due to the dissipation of the mediating field is 
given by 


2. Dissipative Dynamics 


Th e sing l e-mo de h eat baths, characterized by Eqs. 
(42b), (42c), (42e) and (42f), correspond to the unmatched 
side bands of the two pumps, see Eig. These photons 
only carry information about the positions of single oscil¬ 
lators and affect the mechanical modes separately. Pho¬ 
tons in the overlapping side bands partially contribute 
to the coherent exchange of excitations through Jps, but 
also leak out of the cavity, affecting the collective mode 
B. Effective temperatures are defined by the occupa¬ 
tions fii and n. Parameter regimes where or ft become 
negative are to be understood in the sense of statisti¬ 
cal mechanics and indicate anti-damping of the mode in 
question. 

Eor pumps that are arranged symmetrically around the 


F = Fn -h Pini + F 2 n 2 (43) 

and is plotted in Fig. |^as a function of A and n. Each 
contribution is a Lorentzian with width /A centered on 
the blue detuned side of the cavity resonance at A = —Q 
and A = —LI ± 6uj. Eor large the ratio J/F is typi¬ 
cally smaller than unity for A in the vicinity of LI and 
thus one would have to conclude that in this regime high- 
fidelity transfer of quantum states between oscillators is 
not possible. However, for A ^ (2 this ratio becomes a 
monotonically increasing function of A and can in prin¬ 
ciple be made arbitrary large [42]. This remains true 
irrespective of the cavity’s quality factor and is analo¬ 
gous to atomic physics, where the decoherence due to 
the spontaneous emission from an excited state can be 
circumvented by coupling two lower lying states via a far 
detuned two-photon transition. 












we find 


(a) 

noj) 



[DtU=[D^],,= [Dt];,= [D^]l, 

^ -G 1 G 2 

k/2 + i{A - 5 lo/2) 

ra2i = rai2=ra;, = M:, 

^ -G 1 G 2 

Kjj^l -|- -|- 5uj 


( 45 ) 


(46) 


The resonant terms in the matrices Df contribute to 
dissipation 


FIG. 5: (Color online) Sketch of the photon power spectral 
density neglecting the cavity profile, (a) and (b) show the 
power spectral density due to the two pumps in red (solid) 
and blue (dashed) respectively. The two pumps, indicated 
as arrows at ujli and ojl 2 ^ are arrang ed i n the parametric 
amplifier configuration specified by Eq. (28). We emphasized 
the overlapping sidebands captured by Eqs. (44) with dashed 
lines. 


B. Two-mode Parametric Amplification 


Let us now turn to the case where the resonant cou¬ 
pling results in two-mode parametric amplification as de¬ 
scribed by Hamiltonian (§. This interaction leads to fast 
growth of the amplitudes of both modes. The phase of 
each oscillator separately becomes uniformly distributed, 
while the phase of their collective motion becomes in¬ 
creasingly well defined. In the absence of dissipation, the 
resulting two-mode state is not separable and exhibits 
squeezing of the uncertainty of one collective quadratures 
of the system below its zero-point value m- Addition¬ 
ally this systems represents a pair of positive and neg¬ 
ative mass oscillators containing a quantum-mechanics 
free subsystem [36]. Even for degenerate modes, the in¬ 
teraction © needs auxiliary fields to become resonant. 
For the condition (28) we easily find 


— (jj 2 — A 2 T — At 5uj j^l (44a) 
Ai — cji = A 2 T ci;2 = A — 5uj/2. (44b) 


A contributing resonant scattering process is given in 
Fig. 0 (c) and the side band configuration in frequency 
space is sketched in panel (b) of Fig. The pumps are 
arranged such that a blue side band of oscillator 1 is 
overlapped with the red side band of oscillator 2 and 
vice versa. Similar “crossing of side bands” was found in 
other set-ups to lead to entanglement [35| , as is typically 
associated with parametric amplification. 

The resonant coupling terms are found in the off-diagonal 
elements of the matrices and . Using Eqs. (44), 


[Dt],2 = [^ 3 ] 21 

[T]21 = [^3li2 


G 1 G 2 K 
T (A T 
GiG2/^ 

/^^/4 T (A — ^bS)^ 


(47) 

(48) 


Writing out the master equation using these expressions 
and separating unitary from dissipative evolution, we ar¬ 
rive at 

^ = - i [Heff, /9n,] + f _£(S)p^ + f+£(Bt)p^ 

2 

+ ((% + l)^(^i)Pm +n£(6j)/9ni) , (49) 

i=i 


where the effective Hamiltonian is 

^eff = JvK (iib2 T ^ 1 ^ 2 ) 5 (50) 

with 

T — T f GiG2(/^T2iA) ^ 

“ Uv4 + (5a;2/4 - A2 + mA j ’ ^ ^ 

and the collective mode is now 



Again, the commutation relation of this mode depends 
on the relative coupling strengths and vanishes for Gi = 
G 2 , indicating a quantum-mechanics free subsystem [36] . 
Since B involves annihilation and creation operators, one 
cannot unambiguously assign an effective temperature 
to its reservoir. Instead, we will discuss the dynamics 
resulting from these dissipative terms from a numerical 
simulation later on. The parameters characterizing the 
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parametric amplification can be found from Figs. 0(c) 
and A photon is transferred from one pump to the 
other by the creation (or annihilation) of a phonon in 
each mechanical mode. This also clarifies the creation 
of entanglement by this interaction, as phonons in the 
two oscillators are created in pairs and the noise affect¬ 
ing each oscillator becomes correlated. Just as in the 
beam splitter configuration there are necessarily two dis¬ 
tinct pathways for the coherent process separated in fre¬ 
quency by Scj. Consequently systems with photon co¬ 
herence times longer than the inverse of the frequency 
difference, k < Suj, exhibit non-zero detunings where the 
pathways interfere destructively and lead to Jpa = 0. 
These detunings are drawn as lines in Fig. 

For non-vanishing Jpa the unitary interaction corre¬ 
lates the motion of both oscillators while increasing their 
energy exponentially. Whether the mediated interaction 
can entangle the two mechanical modes depends on the 
ratio of Jpa to the total decoherence rate in the system. 

2. Dissipative Dynamics 


FIG. 6: (Color online) Normalized parametric driving 
strength as a function of central detuning and side band res¬ 
olution. The line gives parameters with Jpa = 0. Other pa¬ 
rameters are Gi = G 2 — 0.01 D. Note the formal equivalence 
up to a different normalization compared to Fig. 


dissipative evolution are 

- _ AGlK.A{2n^6uj/2) 

^ ~ [f + (A - 20 - fi)2] [f + {A + 2n+ f )2] ’ 

(53a) 

- _ 4GlKA{2n - 5 lo/2) 


[^ + (A - 20 + ^1^)2] [/^ + (A + 2 O - ^)2] ’ 




(53b) 

r± = 

GiG2/^ 

(53c) 

/t2/4+ (A±(5w/2)2’ 
K‘^/A+{A-2n-5uj/2f 

ni = 

(53d) 

4A(2fi + (5w/2) 

n2 = 

K‘^/A+{R-2Q. + 5uj/2f 

(53e) 

4A(2n - 5uj/2) 


In contrast to the beam splitter case, the collective dy¬ 
namics contained in Eqs. ( [sT] ) and (53c) are independent 
of O and depend on the frequency difference 6uj only. 


1. Unitary Evolution 

The parameter Jpa is proportional to the gain of the 
parametric amplifier and formally equivalent to Jps with 
the substitution Ct ^ 5uj/2. It reaches its maximum Jmax 
at A = + n) and its qualitative behavior is plot¬ 

ted in Fig. The single-quanta processes leading to 


For the parametric amplifier the dissipative behavior 
is somewhat less intuitive than in the beam splitter situ¬ 
ation. While the effect of the two pseudo-thermal, single¬ 
mode reservoirs leads to decoherence of each oscillator’s 
quantum state with rate Vi{ni + 1 ) towards a thermal 
state with occupations the collective mode B is cou¬ 
pled to a non-trivial reservoir. To understand its phe¬ 
nomenology we solved the master equation 

^ = V+C{B)p^ (54) 


numerically for the initial vacuum state. The probability 
distributions for the two single mode oscillators and one 
of their collective modes after t = (Fa)”^ are depicted in 
Figs, [^together with half-maximum contours compared 
to their initial state. The dissipation term (54) increases 
the variance of both quadratures for each oscillator sep¬ 
arately, albeit at a different rate. On the collective mode 
it heats only a single quadrature, leaving the orthogo¬ 
nal one unchanged. These dynamics are the result of a 
back-action evading measurement on the collective mode 
B. For A 7 ^ 0, this evolution is accompanied by the 
diffusion of the orthogonal quadrature due to the non¬ 
vanishing value of r_ in addition to the unitary evolution 
of parametric amplification. 

Just as in the beam splitter case, a cavity with small 
dissipation can resolve the different dissipation channels. 
If the cavity decay rate is smaller than the sum of the 
mechanical frequencies, k < 2Q, one can choose whether 
the s ystem coup les to the single-mode reservoirs through 
( 53a ) and or to the collective, non-trivial reservoir 

(53c). If the decay rate is furthermore smaller than the 
frequency difference, n < 5 either one of the two dif¬ 
ferent collective dynamics governed by Eq. (53c) can be 
chosen separately. 
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FIG. 7: (Color online) Phase-space probability distributions (Wigner functions) for: (a): oscillator 1, (b): oscillator 2, (c): 
collective mode under evolution of given by Eq.(54). The dashed black line is the contour of half maximum, while the thick 
white one marks the same contour for a vacuum state. 



FIG. 8: (Color online) Effective coherences ^pa (as dehned 
in Eq. (55)) as a function of central detuning A and cavity 
decay rate /c.(a) and (c) are for widely separated oscillators, 
6uj/Ct = 1.9, while (b) and (d) are for 6uj/Ct = 0.1. The 
red line denotes the value ^pa = 1. Other parameters are 
Gi = G2 = 0.05 Q. 


3. Entanglement 

The orientation of the squeezed collective quadrature 
in phase-space is not aligned with the squeezing from two¬ 
mode parametric amplification of the effective Hamilto¬ 
nian which occurs at an angle of 7 r /4 with respect to 
the axes chosen in Fig. 0 (c)ian]. In general this dissipa¬ 
tive behavior will therefore counteract the entanglement 
generated by the unitary evolution. The relevant quan¬ 
tity which determines if the two oscillators can be driven 
into an inseparable state is the ratio 

^ rini + r2n2 + r+ + r_' 

Two oscillators can be driven into an entangled state for 
ip A > 1- We plotted this quantity in dependence of n 
and A in Figs, [^for different parameter regimes. 

The behavior of ^pa is qualitatively similar to the 


quantum/classical transition of linearly coupled oscilla¬ 
tors found in [23]. Perhaps most importantly, ^pa scales 
linearly with the central detuning A for high optical de¬ 
cay rates n and can in principle be made arbitrarily large. 
For very large detunings this requires high pump powers, 
since the effective coupling decreases due to smaller in¬ 
tracavity photon numbers, G 1 G 2 oc 1 /A^, and weaker 
effective coupling JpA oc 1 /A, see Eq. 

Just as in the beam splitter case the ratio of unitary 
coupling strength to decoherence rate is analogous to 
two-photon transitions in atomic physics m- For large 
cavity-detunings the lossy photonic state becomes only 
virtually excited and thus decoherence from its dissipa¬ 
tion can be made arbitrarily small. This allows to create 
non-classical correlations in a cavity with a low quality 
factor for sufficiently large central detunings and suffi¬ 
ciently high pump-powers. 


The dissipation rate in the denominator of Eq. (55) as¬ 
sumes that information contained in the photons emitted 
by the cavity is discarded. Using the outcomes of mea¬ 
surements performed on these photons to post-select data 
- or alternatively using these photons directly to control 
a feed-back mechanism - would allow to create entangled 
mechanical modes even in regimes where ip a < 1 . A de¬ 
tailed study of such mechanisms is beyond the scope of 
this work. 


C. Thermal Heat bath 


Throughout this article we have neglected the inde¬ 
pendent thermal heat baths for each oscillator. The 
derivation of the master-equation relies on the coupling 
strength of each oscillator to its bath 7 ^ being much 
smaller than the cavity decay rate k; a condition typi¬ 
cally satisfied for optomechanical systems [1]. The dis¬ 
cussed coherent effects are negligible if thermal decoher- 
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ence rates dominate the system, 


Ui^i > 


G1G2 

K 


(56a) 


where Ui denotes the occupation of oscillator iths ther¬ 
mal reservoir. Even for such thermal systems, however, 
the resonant coupling can significantly alter the dynam¬ 
ics if the effective coupling is at least of the order of the 
reservoir couplings. 


7i - 


G1G2 

K 


(56b) 


These conditions are not necessarily restrictive because 
each Gi can be tuned via the intracavity photon numbers 


IV. SUMMARY 

We investigated a scheme to couple non-degenerate 
mechanical modes resonantly using optomechanical cou¬ 
pling to a single cavity mode. Applying a bi-chromatic 
pump leads to time-dependent optical springs that can 
tune couplings between non-degenerate modes into res¬ 
onance. We derived the reduced master equation for 
the two oscillators containing frequency shifts, unitary 
coupling and non-unitary effects such as dissipation and 
back-action. It is generally applicable for optomechan¬ 
ical systems in the weak coupling regime, ranging from 


macroscopic oscillators to ultracold atomic systems. We 
applied the master equation to discuss two resonant cases 
in depth: a beam splitter and a parametric amplifier. 
While the former allows more general cooling schemes 
and quantum state transfer between different frequency 
regimes, the latter could realize entanglement of spatially 
separated mechanical modes. We showed that these two 
fundamental quantum operations can be executed coher¬ 
ently irrespective of the cavity quality factor. Addition¬ 
ally, both cases exhibit mechanical quantum interference 
for sufficiently low cavity decay rates. 

The coupling scheme allows customized interactions 
among distinct elements in arrays consisting of more than 
two mechanical modes m and the derived master equa¬ 
tion is easily generalized to include more modes. Be¬ 
yond mechanical excitations, the presented theory can 
be adapted to describe a variety of hybrid quantum sys¬ 
tems, for instance ensembles of spins interacting with a 
single cavity mode [38] or spatially separated nitrogen 
vacancy centers in diamond coupled by a high-frequency 
phonon mode [39] . 
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Appendix A: Nonlinear corrections to the resonance 
conditions 

Adopting the notation of Sec. [ll| with uji and 5VLi de¬ 
noting the bare mechanical frequencies and their optical 


spring shift respectively, the two modified resonance con¬ 
ditions read 


Al — A2 — Suj T — 6^2 

Al - A 2 =2Q -b SQi + 6^2^ • (A2) 

These resonance conditions are nonlinear equations with 
respect to the two pump frequencies because the indi¬ 
vidual springs 6fti depend on the position of each pump 
with respect to the cavity resonance. While the individ¬ 
ual shifts are subtracted from each other in the case of 
a beam splitter coupling, they add up for the parametric 
amplifier. For mechanical modes with identical couplings 
they remain small corrections to the bare frequency dif¬ 
ference/sum, but since they provide limits to the fidelity 
of coherent state transfer or the gain of parametric ampli¬ 
fication and the resulting entanglement, these effects have 
to be taken into account for applications in the quantum 
regime. The extent of the correction is measured by the 
dimensionless quantities 


SQi - 6^2 
6uj 

(5Qi T ^^2 
20 ’ 


(A3) 

(A4) 


for the beam splitter and parametric amplifier respec¬ 
tively. We have calculated these expressions numerically 
and plotted them as a function of the bare mechanical 
frequency difference 6uj and the average of the pump de¬ 
tunings in Fig. 
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FIG. 9: (Color online) Corrections of the resonant pump detunings from their bare value. Top row for the beam splitter setting 
, bottom row for parametric amplification, All plots have Qjak = O.OIQ and the optical dissipation rate is 

tv — lOfi for (a) and (d), k. — Cl for (b) and (e) and tv = O.lCl for (c) and (f). 

























